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ABSTRACT 

We study the clustering properties of groups and of galaxies in groups in the DEEP2 Galaxy Redshift 
Survey dataset at z ~ 1 in three separate fields covering a total of 2 degrees 2 . Four measures of two- 
point clustering in the DEEP2 data are presented: 1) the group correlation function for 460 groups 
with estimated velocity dispersions of a > 200 km s _1 , 2) the galaxy correlation for the full DEEP2 
galaxy sample, using a flux-limited sample of 9800 objects between 0.7 < z < 1.0, 3) the galaxy 
correlation for galaxies in groups or in the field, and 4) the group-galaxy cross-correlation function. 
Our results are compared with mock group and galaxy catalogs produced from ACDM simulations. 
Using the observed number density and clustering amplitude of the DEEP2 groups, the estimated 
minimum group dark matter halo mass is M m i n ~ 6 x 10 12 /i _1 A/ Q for a flat ACDM cosmology 
with as = 0.9. Groups are more clustered than galaxies in the DEEP2 data, with a relative bias of 
b = 1.17 ± 0.04 on scales r p = 0.5 — 15 h^ 1 Mpc. Galaxies in groups are also more clustered than the 
full galaxy sample, with a scale-dependent relative bias which falls from b ~ 2.5 ± 0.3 at r p = 0.1 h^ 1 
Mpc to b ~ 1 ± 0.5 at r v = 10 Mpc. The correlation function of galaxies in groups has a steeper 
slope (7 ~ 2.12 ± 0.06) than for the full galaxy sample (7 ~ 1.74 ± 0.03), and both samples can be fit 
by a power-law on scales r p — 0.05 — 20 h^ 1 Mpc. We empirically measure the contribution to the 
projected correlation function, w p (r p ), for galaxies in groups from a 'one-halo' term and a 'two-halo' 
term by counting pairs of galaxies in the same or in different groups. The projected cross-correlation 
between group centers and the full galaxy sample, which is sensitive to the radial distribution of 
galaxies in and around groups, shows that red galaxies are more centrally concentrated in groups 
than blue galaxies at z ~ 1. DEEP2 galaxies in groups appear to have a shallower radial distribution 
than that of mock galaxy catalogs made from N-body simulations, which assume a central galaxy 
surrounded by satellite galaxies with an NFW profile. Using simulations with different halo model 
parameters, we show that the clustering of galaxies in groups can be used to place tighter constraints 
on the halo model than can be gained from using just the usual galaxy correlation function alone. 
Subject headings: galaxies: statistics — cosmology: large-scale structure of universe 



1. INTRODUCTION 

Groups of galaxies populate an intermediate range in 
density-contrast between galaxies and clusters and oc- 
cupy a regime that is critical to understanding hier- 
archical galaxy formation in ACDM models. Merger 
events between galaxies likely occur within groups rather 
than clusters due to their lowe r velocity dispersions (e.g., 
iBarnesI ll98H: lOstrikerl Il980l) . Galaxy groups should 
also be more easily related to dark matter halos than 
galaxies themselves, which can have a complicated halo- 
occupation function that depends significantly on galaxy 
properties. In order to understand and test galaxy for- 
mation and evolution models it is useful to relate galax- 
ies to observable groups as a proxy for their parent dark 
matter halos. The cu r rent halo model paradigm (e.g., 
i Coorav & Shethl 12001 jKravtsov et aT1l200l iMa fc Frvl 
l2000UPeacock fc Smithl2000l:lSeliakl2000(l provides a sta- 
tistical analytic measure for relating galaxies to their 
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dark matter halos. A key statistic in the halo model 
is the halo occupation distribution (HOD) which mea- 
sures the probability of a halo of a given mass hosting 
N galaxies. The halo model also naturally explains the 
small deviations seen in the clustering of galaxies at z ~ 
from a power-law model, where there is a transition from 
galaxies within a single halo and between different halos. 
The clustering of groups and group galaxies depends not 
only on halo model parameters but also on the nature of 
bias and the details of hierarchical structure formation, 
as well as cosmological parameters such as Q, m and as ■ 

Observationally, clusters of galaxies have b een shown 
to be very strongly clustered l|Bahcalll IT988) . with the 
clus tering strength depending on the richness of the clus- 
ter. iKaiserl l)1984f) show that the large clustering scale- 
length of massive and rare Abell clusters can be explained 
by a simple model in which these clusters formed in re- 
gions where the primordial density enhancement was un- 
usually high. Objects forming in the densest peaks would 
naturally be biased tracers of the underlying dark matter 
field such that massive clusters would have a higher corre- 
lation amplitude than that of galaxies. Since their masses 
are intermediate, galaxy groups are therefore expected to 
have clustering properties between those of galaxies and 
clusters. 

The first papers analyzing the clustering of groups 
in local redshift surveys at z ~ present con- 
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flicting results and were hamp ered by small sam- 
ples and cosmic variance fe.g.. Uing fc Zhanel Il988t 
i Maia fe da Costal fl990t iRamella. Geller. fc Huchral 
^)9oTh ^^Isarti-Battistoni. Invernizzi^^Bonomettol 
( 1997) investigated the effect of changing the linking- 
length parameters in the Friends-of-Friends (FoF) 
algorithm and its effect on the clustering signal 
and found that these early papers used too large 
a linking length, which led to a diminished clus- 
tering strength due to the presence of i nterl opers . 
iTrasarti-Battistoni. Invernizzi. fc Bonomettol l)1997f) 
found in their data from the Perseus-Pisces redshift 
survey, using two fields and ~ 50 and 200 groups in each 
field, that groups are approximately twice as clustered 
as galaxies, but with significant error bars. All of these 
papers showed that determining the clustering proper- 
ties of groups is a tricky endeavor, which can depend 
quite sensitively on the volume and magnitude depth of 
the survey, the handling of the selection function and 
varying completeness, and the method used to identify 
groups. In addition, these analyses all suffered from 
significant uncertainties, from both Poisson statistics, 
due to the small number of groups in the surveys, and 
cosmic variance, due to the small volumes surveyed, 
which was not quantified in any of these papers. 

Significant advances h ave recently been made with 
much larger datasets. iGirardi. Boschin. fc da Costal 
(2000), with the combined CfA2 and SSRS2 surveys, 
have a sample of 885 groups in a volume of ~ 3 x 
10 5 /i~ 3 Mpc 3 , much larger than earlier surveys. With 
this large sample size, they are able to construct volume- 
limited subsamples and investigate the dependence of 
clustering on group properties, finding that groups with 
more members and/or larger internal velocity dispersion 
are m ore strongly clustered. iMerchan. Maia. fc Lambasl 
(2000) use a sample of 517 groups from the Updated 
Zwicky Catalog and 104 groups from the SSRS2 to show 
that groups are at least twice as clustered as galaxies, and 
that more massive groups have larger clustering strength. 

Most recently, the 2dF Galaxy Redshift Survey has 
provided a vast dataset with which to study large-scale 
struct ure at z < 0-2. Using data fr o m the 100k re- 
lease, iZandivarez. Merchan. fc Padillal l)2003(l measure 
the clustering of groups as a function of virial mass and 
find that more massive groups are more clustered and 
that their measurements match the clustering of dark 
matter halos in a ACDM N-body simulation . Using 
the completed 2dF survey, iPadilla et al.l l)2004f ) analyze 
group clustering as a function of luminosity and show 
that while the least luminous groups actually cluster less 
than the galaxies in the survey, there is a strong relation 
between group luminosity and correlation length, and 
the most luminous groups (with L~4x 10 11 /i _2 L Q ) are 
~ 10 times more clustered as the least luminous groups 
(with L ~ 2 x 10 10 /i~ 2 L Q ). The relation between clus- 
tering scale length and mean group separation that they 
find continues the trend seen on larger scales for clus- 
ters. They find very good agreement between their data 
and mock catalogs constructed from ACDM simulations 
and semi-analytic galaxy evolution recipes. T hese same 
conclusions are reached bv lYang et alJ l)2005c|) . who also 
use the 2dF data to measure the clustering of groups as 
a function of luminosity. It appears that at z ~ there 



is now convergence among group clustering analyses. 

The extensive 2dF group catalogs have now also al- 
lowed studies of groups beyond simple measures of the 
correlation function of groups. Several authors measure 
the radial profile of gala x ies in group s usin g 2dF data 
( Colli ster fc Lahavl 120051: iDiaz et all 120051: lYang et alJ 
2005b), and analyze t he cross-correlation between group 
centers and galaxies llYang et al.ll2005b|) . These papers 
find that galaxies in groups are less concentrated than 
dark matter particles in simulations, and also find that, 
locally, the centers of groups are preferentially populated 
by red galaxies compared to blue galaxies. The HOD has 
now also been measured directly at z ~ using counts of 
galax i es in groups of diffe rent masses (jCollister fc Lahavl 
120051 lYang et aUl2005a|) . constraining halo model pa- 
rameters locally. Clustering measures of the correlation 
function of all galaxies (not just those in groups) at inter- 
mediate redshifts indicate that the H OD does not change 
significantly between z ~ 1 and z ~ (Ph leos et al . 2005; 
lYan. Madgwick. fc Whitell200^ . 

These measurements have only been performed with 
local samples; group catalogs have not been available at 
intermediate- or high-redshift. In this paper we present 
the first analysis of group clustering at z ~ 1, using 
group catalogs from the DEEP2 Galaxy Redshift Survey 
(|Oerke et all 120051) . The high resolution of the DEEP2 
data allows us to identify groups in three dimensions re- 
gardless of their galaxy properties, using only their over- 
density in space. We focus here on the clustering of 
groups and galaxies within groups at z ~ 1 , as these mea- 
sures, when combined with similar measures at z ~ 0, 
will provide constraints on galaxy evolution and struc- 
ture formation models. We also show how these measures 
can constrain the HOD at z ~ 1. 

Four measures of two-point clustering in the DEEP2 
dataset between 0.7 < z < 1.0 are analyzed: 1) the group 
correlation function, 2) the galaxy correlation function 
for DEEP2 galaxies, 3) the galaxy correlation function 
for galaxies in groups, and 4) the group-galaxy cross- 
correlation function. These clustering measures can be 
used to constrain the halo model parameters by com- 
paring the data to mock catalogs with different HODs. 
The first clustering measure, when combined with the ob- 
served number density of groups in our sample, is used 
to estimate the typical dark matter masses of the halos 
the groups studied reside in. The second measure pro- 
vides constraints on the HOD, though not in an entirely 
unique way; further constraints on the HOD are provided 
by the last two measures. For the clustering of galaxies 
in groups, we empirically distinguish between 'one-halo' 
and 'two-halo' terms using pairs of galaxies within the 
same group and pairs in different groups, respectively. 
The clustering of galaxies in groups and the group-galaxy 
cross-correlation function can also constrain the radial 
distribution of galaxies within groups when compared 
with mock catalogs. 

Lastly, groups may be used to constrain cosmological 
parameters like the dark energy equation of state, w, if 
their bivariate distribution in redshift z a nd velocity dis- 
persi on a can be accurately measured ((Newman et alJ 
2002). This test, however, requires that the relation be- 
tween group velocity dispersion and dark matter halo 
mass be known a nd accurately calib r ated. It has recently 
been suggested ijLima fc Hul 120041: iMaiumdar fc Mohrl 
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120041) that the clustering properties of galaxy clusters 
may be used for "self-calibration", since the clustering 
properties of halos can be predicted as a function of mass 
and compared to the measured clustering. The group 
correlation function results presented here should be use- 
ful as such a self-calibration procedure for future DEEP2 
studies. 

An outline of the paper is as follows: §2 briefly de- 
scribes the DEEP2 Galaxy Redshift Survey and the sam- 
ple of galaxy groups used here, as well as the mock 
galaxy catalogs constructed for the survey. §3 discusses 
the methods used to calculate the two-point correlation 
functions. We present clustering results for groups in 
§4, where we compare with simulations and estimate the 
minimum dark matter halo mass for our groups. In §5 
we analyze the clustering of the full galaxy sample and 
galaxies in groups in the DEEP2 data and in mock cata- 
logs and show the contribution to the correlation function 
for galaxies in groups from the 'one-halo' and 'two-halo' 
terms. §6 presents the cross-correlation between the full 
galaxy sample and group centers, which depends upon 
the radial profile of galaxies in groups. The relative bi- 
ases between galaxies in groups and all galaxies and be- 
tween groups and galaxies are presented in §7. Mock 
catalogs with different halo models are used to illustrate 
how these various clustering measures can be used to con- 
strain parameters of the halo model in §8. We conclude 
in §9. 

2. DATA SAMPLE AND MOCK CATALOGS 

In this paper we use data from the DEEP2 Galaxy 
Redshift Survey, which i s an ongoing proje ct using the 
DEIMOS spectrograph l|Faber et all 12002ft on the 10- 
m Keck II telescope to survey optical galaxies at z ~ 
1 in a comoving volume of approximately 5xl0 6 h~ 3 
Mpc 3 . Using ~ 1 hr exposure times, the survey will 
measure redshifts for ~ 40, 000 galaxies in the red- 
shift range 0.7 ~ z ~ 1.5 to a limiting magni tude of 
i?AB = 24.1 ijCoil et alJl2004^lFaber et al.ll200l . Spec- 
troscopic targets are pre-selected using a color cut in 
B — R - R — I space to ensure that most galaxies lie 
beyond z ~ 0.75. This color-cut results in a sample 
with ^90% of the targeted objects at z > 0.75, miss- 
ing only ~3% of the z > 0.75 galaxi es which meet 
our magnitude limit l)Davis et alJ [20021. Due to the 
high dispersion (R ~ 5, 000) of our spectra, our red- 
shift errors, determined from repeated observations, are 
~ 30km s _1 . Restfr ame (U — B)n colors have been de- 
rived as described in l)Willmer et alJ2 005). Details of the 
observations, ca talog cons t ruction and data reductio n 
can be found inlCoil et al I ll2004bl): iDavis et all <2Q02fl : 
Davis. Gerke. fc Newman! l)2004|) ; iFaber et alJ l|2005|) . 

The completed survey will cover 3 square degrees of 
the sky over four widely separated fields to limit the im- 
pact of cosmic variance. Each field is comprised of two 
to four contiguous photometric 'pointings' of size 0.5 by 
0.67 degrees. Here we use data from six of our most 
complete pointings to date, in three separate fields. We 
use data from pointings 1 and 2 in the DEEP2 fields 
2, 3 and 4. Each DEEP2 pointing corresponds to a 
volume of comoving dimensions ~ 20 x 27 x 550 hr 1 
Mpc in a ACDM model for 0.7 < z < 1.0. The 
total volume of the sample used in this paper is thus 
~ 1.8 x 10 6 /i~ 3 Mpc 3 . To convert measured redshifts 



to comoving distances along the line of sight we assume 
a flat ACDM cosmology with Q m — 0.3 and J7a = 0.7. 
We define h = i7o/(100 km s" 1 Mpc" 1 ) and quote cor- 
relation lengths, ro, in comoving h^ 1 Mpc. Throughout 
the paper, we quote empirical errors calculated from the 
variance across the six pointings. In so doing we have 
treated the six DEEP2 pointings as being entirely inde- 
pendent, and they are not: there are two adjacent point- 
ings in each of three independent fields. We estimate 
from Monte Carlo simulatio ns using the methods out- 
lined in lNewman et all l)2002|) that, for the distribution of 
pointings used here, the measured standard error should 
be increased by 17 ± 5%, for errors which are dominated 
by cosmic variance; this is a conservative assumption, as 
(especially for the group-group correlations) other con- 
tributions such as shot noise are significant. 

A description of the methods empl oyed to detect 
groups is presented in lGerke et alJ l)2005[K along with de- 
tails of the group cata log. A Voronoi-Delau nay Method 
(VDM) group-finder l|Marinoni et all I2002T) is used to 
identify galaxy groups. This method searches for galaxy 
overdensities in redshift space, using an asymmetric 
search window to account for redshift-space distortions. 
The advantage of the method over traditional FoF group- 
finding methods lies in the fact that it has no fixed length 
scale, but instead uses an adaptive search radius based on 
estimated group richness. The VDM group-finder thus 
avoids a common problem of FoF methods, in which clus- 
tered groups are merged together along filaments. 

We use three galaxy samples in this paper; all are 
drawn from the same volume as the group sample. The 
full galaxy sample includes a total of 9787 galaxies in 
the redshift range 0.7 < z < 1.0 in the same six point- 
ings used for the group analysis. We also split the full 
galaxy sample into subsamples of field and group galax- 
ies, with 5947 and 3840 objects in each. Throughout this 
paper we use the terms "all galaxies" and "the full galaxy 
sample" interchangeably; they include both the field and 
group galaxy samples. 

Our main group sample consists of all DEEP2 groups 
with an estimated velocity dispersion a > 200 km s , 
with a total of 460 groups and ~ 50 — 100 groups per 
pointing. We do not make a distinction between groups 
and clusters - clusters are simply the larger groups; they 
are included in all analyses, but have minimal impact 
on this work due to their rarity. The a > 200 km s _1 
group sample has similar completeness (74 ± 5%) and 
purity (57 ± 3%) as samples with higher a cutoffs, as 
estimated using mock catalogs described below. Here 
completeness is defined as the fraction of real groups (de- 
fined as a group of galaxies belonging to the same parent 
halo, where halos are identified using FoF in real space 
on the dark matter particles in the simulations) that 
are successfully identified in the recovered group catalog 
created by the group-finder, and purity is the fraction 
of recovered group s that correspond to real groups (see 
IGerke ei ^ al.| l|2Q Q5l) for more details). As shown in Fig. 
8 of IGerke et alJ 1)20051 ). these statistics are nearly in- 
dependent of group velocity dispersion. There are more 
recovered groups than real groups in these mock catalogs 
by a factor of 1.4. The a > 200 km s _1 group sample has 
a relatively high "galaxy success rate" of 70 ±1%, defined 
as the fraction of galaxies in real groups that are iden- 
tified as group galaxies in our recovered group catalog, 
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and an interloper fraction of 43 ± 1 %. We do not use th e 
a > 350 km s _1 cut that was used in lOerke et alJ |2005), 
as that cut was used to define a sample that recovered 
different physical properties than the ones relevant here; 
instead of the distribution of groups in redshift and ve- 
locity dispersion, the relevant parameters here are the 
positions of groups and the interloper fraction (the frac- 
tion of identified group galaxies that are not actually in 
groups). 

The observed richness distribution of our groups in the 
DEEP2 data is roughly a power-law, with most groups 
having two observed galaxies; the largest groups have 
~ 10 galaxies, though there are only a few groups this 
large. Fig. ^ shows the observed richness distribution 
for groups in the DEEP2 data as well as for real and 
recovered groups in our mock catalogs, for groups with 
0.7 < z < 1.0 and a > 200 km s" 1 . 

A histogram of the redshift distribution of our group 
sample is shown in Fig. [21 We restrict the analyses here 
to the redshift range 0.7 < z < 1.0 to minimize sys- 
tematic effects. While groups are found in the survey 
to higher redshifts, those groups are likely to be more 
massive and hence not as representative of groups in our 
sample as a whole. Additionally, the i?-band target se- 
lection of the survey corresponds to a bluer restframe 
color-selection at higher redshift; this results in fewer red 
galaxies bein g targeted at higher redshift compared to 
blue galaxies l)Willmer et al.ll2Q05^l . which could system- 
atically affect the group richness and velocity dispersion 
estimates at z > 1.0. The spatial distribution of DEEP2 
groups in three of the six pointings used is shown in Fig. 
[3J Shown are pointings 2 in the DEEP2 fields 2, 3 and 
4. 

The mock galax y catalogs used throughou t the pa- 
per are described in lYan. White, fc Coill (|2004|) : relevant 
details are repeated here. The mock catalogs are con- 
structed from N-body simulations of 512 3 dark matter 
particles with a particle mass of m = 1 x 10 10 h~ 1 M© 
in a box with dimensions 256 /i -1 Mpc on a side, for a 
ACDM cosmology with cr 8 = 0.9. Dark matter halos 
were identified using FoF in real space and galaxies were 
placed within halos using a halo model prescription. To 
populate dark matter halos with galaxies, two functions 
need to be specified. The first is the halo occupation dis- 
tribution function (HOD), which is the probability that 
a halo of mass M hosts N galaxies, P(N\M). The second 
is the spatial distribution of galaxies within halos. The 
first moment of the HOD function, the average number 
of galaxie s as a function of halo m ass M, is shown in 
Fig. 1 in lYan. White, fc Coill ^2004) and in Fig. 12 of 
this paper (labeled as "B256") for the HOD used in the 
mock galaxy catalogs here. This HOD is constrained at 
z ~ with the 2dF luminosity function and luminosity 
dependent two-point correlation function of galaxies and 
at z ~ 1 with the DEEP2 £(r) and the COMBO-17 lu- 
minosity functions. In §8 of this paper we compare clus- 
tering results for mock catalogs with two different halo 
models; for the bulk of the paper, h owever, we use light- 
cones with the HOD as published in lYan. White. &: Coill 
(2004). Once the number of galaxies in each halo and 
their corresponding luminosities are known, the most lu- 
minous galaxy is assigned to the center of mass of the 
halo, and positions and velocities for the other galax- 
ies are drawn randomly from those of the dark matter 



particles. No other radial or velocity bias is included in 
assigning galaxies to particles, so that galaxies trace the 
mass and velocity distributions of the dark matter par- 
ticles in the halo. The spatial distribution of galaxies 
follow the dark matter density profile, which on average 
is an NFW profile. Although no velocity bias is included, 
the velocity dispersions of the galaxies are systematically 
smaller than those of the mass (see Fig. 8 in Yan, White 
& Coil 2004), due to the fact that the most luminous 
galaxy is always assigned to the particle at the center 
of mass. Here we use a set of twelve independent cata- 
logs which each have the same spatial extent as a single 
DEEP2 pointing (~ 0.5 by 0.67 degrees). 

The center of each DEEP 2 group (needed for the clus- 
tering measures performed here) is measured as the me- 
dian of the positions (in comoving x, y, and z) of the 
galaxies identified in that group. Errors on the positions 
of the recovered groups are estimated using the differ- 
ences in the mock catalogs between the centers of real 
groups and recovered groups. As stated above, there is a 
factor of 1.4 more recovered groups in the mock catalogs 
than real groups for a > 200 km s . In the mock cat- 
alogs 74% of recovered groups have a real group within 
r p = 1 h^ 1 Mpc, while 48% have a real group within 
r p = 0.2 h^ 1 Mpc and 41% have a real group within 
r p = 0.1 h^ 1 Mpc. Here r p is the projected distance 
on the plane of the sky, which is the relevant distance 
for the projected clustering used in §5 and later in the 
paper; in §4 we use the redshift-space correlation func- 
tion, but only measure it on scales s > 2 h^ 1 Mpc. We 
therefore believe that on scales larger than 1 h^ 1 Mpc 
our results are robust to errors in the group positions, 
while on smaller scales, where we will measure the cross- 
correlation between group centers and galaxies, there is 
likely to be some degradation in the signal due to posi- 
tion errors; this is discussed more in §9. However, the 
mock catalogs have been treated in an identical manner 
as the data so that comparisons between the data and 
mock catalogs are unaffected. 

3. METHODS 

The two-point correlation function £(r) is defined as a 
measure of the excess probability above Poisson of find- 
ing an object in a volume element dV at a separation r 
from another randomly chosen object, 

dP = n[l + f(r)]dV, (1) 

where n is the mean n umber density of the object in 
question llPeebleslll980l 

Measuring £(r) requires constructing a random catalog 
with the same selection criteria and observational effects 
as the data, to serve as an unclustered distribution with 
which to compare. For each data sample we create a 
random catalog with the same overall sky coverage and 
redshift distribution as the data. This is achieved by 
first applying the two-dimensional window function of 
our data in the plane of the sky to the random catalog. 
Our overall redshift success rate is gtrsim70% and is not 
entirely uniform across the survey; some slitmasks are 
observed under better conditions than others and there- 
fore yield a slightly higher completeness. This spatially- 
varying redshift success rate is taken into account in the 
spatial window function which is applied to both the ran- 
dom catalog and the mock catalogs, such that regions of 
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the sky with a higher completeness have a correspond- 
ingly higher number of random points or more objects in 
the mock galaxy catalogs. This ensures that there is no 
bias introduced when computing correlation statistics or 
comparing the data to the mock catalogs. We also mask 
the regions of the random and mock catalogs where the 
photometric data have saturated stars and CCD defects. 

We then apply a selection function, (f>(z), defined as 
the probability of observing a group as a function of red- 
shift, so that the random catalog has the same overall 
rcdshift distribution as the data. The selection func- 
tion for groups in the DEEP2 survey is determined by 
smoothing the observed redshift distribution of groups in 
the catalog used here, as shown in Fig. |21 The appar- 
ent overdensity at z ~ 0.85 seen in this figure is due to 
structures in several pointings and does not significantly 
affect our results, which are averaged over all pointings 
and a wide redshift range. As the group catalog includes 
three separate fields, the data are combined when calcu- 
lating <f>(z), which reduces the effects of cosmic variance. 
Using a smoothed redshift distribution to estimate the 
selection function can cause a systematic bias, but we 
estimate that this bias is less than the errors on £ (r) . In 
what follows, for both the data and the mock catalogs 
the redshift range over which we compute £(r) is limited 
to 0.7 < z < 1.0, the redshift range over which the se- 
lection function varies by less than a factor of two. Each 
of the six pointings has between ~ 50 — 100 groups, and 
each of the twelve independent mock catalogs contains 
~ 100 groups. We use ~ 3, 000 random points in each 
pointing to calculate the correlation function for groups 
and ~ 20, 000 random points to calculate the correlation 
function for galaxies. 

The two-point correlation function is measured for 
both groups and galaxies using the ILandv fc Szalavl 
(1993) estimator, 
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where DD, DR, and RR are pair counts in a given sepa- 
ration range in the data-data, data-random, and random- 
random catalogs, and ud and Hr are the mean number 
densities in the data and random catalogs. This estima- 
tor has b een shown to pe rform as well as the Hamilton es- 
timator l|Hamiltonll993D but is preferred as it is relatively 
insensitive to the size of t he random catalog and han- 
dles edge corrections well l|Kerscher. Szapudi. fc Szalavl 
200(J. As the magnitude limit of our survey results in a 
non- uniform selectio n function, a standard J 3 -weighting 
scheme is applied l|Davis fc Huchral l!982|K which at- 
tempts to weight each volume element equally while min- 
imizing the variance on large scales. As the redshift range 
is limited to 0.7 < z < 1.0 for all analyses in this paper, 
this effect is not large (< 20%). 

Measurements of the cross-correlation between two 
samples (presented in §6) use a symmetrized version of 
Eqn. 2. Each data sample, with pair counts labeled 
D\ and Z?2, has an associated random catalog, with pair 
counts R\ and i?2 , with the same selection function as the 
data. After normalizing each data and random catalog 
by its number density, the cross-correlation is estimated 



RiR: 



[DtLh - DiR 2 - D 2 Rx + R1R2] . (3) 



Redshift-space distortions due to peculiar velocities 
along the line-of-sight will introduce systematic effects 
to the estimate of £(?"). At small separations, random 
motions within a virialized overdensity cause an elonga- 
tion along the line-of-sight ("fingers of God"), while on 
large scales, coherent infall of galaxies into forming struc- 
tures causes an apparent contraction of structure along 
the line-of-sight (the "Kaiser effect"). What is actually 
measured then is £(s), where s is the redshift-space sepa- 
ration between a pair of galaxies. As we will show in the 
next section using mock group catalogs, for galaxy groups 
over the scales relevant here, the systematic effects of 
redshift-space distortions are of order 20%. Mock cata- 
logs are used to correct for redshift-space distortions and 
infer £(r) from £(s)by multiplying the observed £(s)for 
DEEP2 groups by the ratio of £(r)to £(s)for groups in 
the mock catalogs. 

However, redshift-space distortions are more signifi- 
cant when measuring £ for galaxies, as smaller scales are 
probed. We are able to measure £ on small scales for 
the galaxy sample both because of the larger sample size, 
which allows us to stably measure the clustering on small 
scales, and also because there is no exclusion radius, un- 
like with groups, which typically have a radius of R ~ 1 
h^ 1 Mpc. To uncover the real-space clustering proper- 
ties of galaxies we measure £ in two dimensions, both 
perpendicular to a nd along the line of sight. Following 
iFisher et alJ (^994), vi and V2 are the redshift positions 
of a pair of galaxies, s is the redshift-space separation 
(vi — va), and 1 = h(vi + V2) is the mean distance to 
the pair. The separation between the two galaxies across 
(r p ) and along (w) the line of sight are defined as 

«= W (4) 



r p = y/ s • s - 7r 2 . (5) 

In applying the ILandv fc Szalavl l|l993h estimator to 
galaxies, pair counts are computed over a two- 
dimensional grid of separations to estimate £(r p , tt). 

To recover £(r) £(r p , tt) is projected along the r p axis. 
As redshift-space distortions affect only the line-of-sight 
component of £(r p ,7r), integrating over the tt direc- 
tion leads to a statistic w p (r p ), whic h is independent 
of red shift-space distortions. Following IDavis fc Peebles! 

GUI, 

poo pOO 

w p (r p ) = 2 dir £(rv,7r) = 2 / dy £(r 2 p + y 2 ) 1 ' 2 , 



(6) 

where y is the real-space separation along the line of 
sight. If is modeled as a power-law, £(r) = (r/ro) -7 , 
then ro and 7 can be readily extracted from the projected 
correlation function, w p (r p ), using an analytic solution to 
Equation [|J] 



w p (r p ) = r p 



r(!)W) 



(7) 



where T is the gamma function. A power-law fit to w p (r p ) 
will then recover ro and 7 for the real-space correlation 
function, £(r). 
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4. GROUP CLUSTERING RESULTS 

This section presents results on the clustering of groups 
at z *~ 1. The two-point correlation function for groups 
in the DEEP2 data is measured. The clustering proper- 
ties of groups in our mock catalogs are analyzed, which 
allows us to quantify our systematic errors and to correct 
the observed clustering in the DEEP2 data for redshift- 
space distortions and the effects of the group-finder. The 
number density of groups in the DEEP2 data is then used 
to infer the minimum dark matter halo mass which hosts 
our groups; we then compare the clustering of halos with 
this mass to the clustering of our group sample. 

4.1. Clustering of Groups in DEEP2 Data and Mock 

Catalogs 

The left panel of Fig. 0] shows the observed two-point 
correlation function in redshift space, £(s), for our a > 
200 km s -1 group sample, for both the full sample (solid 
line) with two or more galaxies in each group (N > 2) 
and for a subsample with four or more galaxies in each 
group (N > 4). Also shown is £(s) for the recovered 
group sample in the mock catalogs with a > 200 and 
N > 2 (dotted line). We show the standard error across 
the six pointings; this empirical error therefore includes 
both cosmic variance and Poisson error. The N > 4 
data sample shows a stronger clustering amplitude than 
the N > 2 data sample. This is likely due to the N > 4 
sample containing more massive groups, on average, than 
the N > 2 sample. The N > 4 sample is also less likely 
to have interlopers, which may also increase the observed 
clustering. 

Power-law fits to £(s) over the range s = 3 — 20 /i -1 
Mpc are given in Table 1. Below s ~ 3 /i -1 Mpc £(s) 
does not continue to rise as a single power-law in our 
data; this lack of pairs on scales s < 3 hr 1 Mpc is likely 
due to the finite physical extent of groups. To take into 
account the covariance between so and 7, we perform 
a x 2 minimization and marginalize over each parame- 
ter separately. This procedure leads to errors which are 
roughly a factor of two larger than if we neglected this 
covariance. 

The mock galaxy catalogs described in §2 are used to 
quantify systematic errors in our group clustering anal- 
ysis. Effects that may be introduced by the group- 
finder are tested by comparing £(s) in redshift space 
for both real and recovered groups. The results are 
shown in Table 1. The recovered groups have a slightly 
higher clustering amplitude (5%) than the real groups. 
Redshift-space distortions are quantified by measuring 
£(s) in redshift space and £(r) in real space for real 
groups. Redshift-space distortions appear to enhance 
the clustering properties of groups by ~ 20% on scales 
s ~ 2 — 15 ft -1 Mpc and decrease the clustering am- 
plitude on smaller scales. The clustering scale length 
increases by ~ 20% when measured in redshift space, 
while there is no change to the slope over the scales used 
here. The effects of our slitmask target selection algo- 
rithm on the clustering of groups are also tested using 
the mock catalogs. Our target selection code determines 
which galaxies would be observed on slitmasks; because 
spectra from neighboring galaxies can not overlap on the 
CCD, not all galaxies can be selected to be observed. In 
particular, in overdense regions on the plane of the sky 



the number of galaxies which can be observed decreases 
in a non-trivial way. By running our slitmask target se- 
lection code on the mock catalogs we can quantify the 
effect on the measured £(r). 

Under the assumption that corrections for the above 
effects should be proportional to the clustering strength, 
these results from the mock catalogs can now be used to 
correct the observed £(s) for groups in the DEEP2 data 
for target selection effects, redshift space distortions, and 
the group-finder, in order to estimate £(r) in real space 
for real groups. We can either apply corrections to the 
power-law fits themselves or to the data points as a func- 
tion of scale. If we apply the corrections to so and 7 as 
measured for the DEEP2 groups, the corrections would 
infer that r Q = 6.8 ± 0.6 hr 1 Mpc and 7 = 1.4 ± 0.2 
for the real-space correlation function of real groups in 
the DEEP2 data with N > 2. If we explicitly correct 
the observed £(s) for the DEEP2 groups as a function of 
scale, using the ratio of £(v)/£( s ) in the mock catalogs, 
the resulting £(r) has a power-law fit of ro = 6.2 ± 0.4 
hr 1 Mpc and 7 = 1.5 ± 0.2, within the la error of the 
inferred values. The corrected £(r) is shown as the solid 
line in the right panel of Fig. 01 for both the N > 2 
and N > 4 samples. The clustering of groups in the 
DEEP2 data is therefore l-3<7 lower than the clustering 
of real groups before target selection in the mock cata- 
logs, where r = 7.4 ± 0.4 and 7 = 1.6 ± 0.2, for N > 2. 

4.2. Minimum Group Mass Derived from Clustering 

Results 

The minimum dark matter halo mass that our groups 
reside in can be inferred from the observed number den- 
sity and clustering strength of our groups using cither 
analytic formulations of the dark matter halo mass func- 
tion or by comparing directly to simulations. Here we 
estimate a minimum dark matter halo mass from the ob- 
served number density using analytic theory and then 
compare the corresponding predicted clustering strength 
of those halos in both theory and simulations with the 
observed clustering strength of our groups. 

The a > 200 km s -1 group sample has an observed 
density of n = 4.5 x 10 -4 ft, 3 Mpc -3 , calculated from 
the observed number of groups between z — 0.75 — 0.85, 
where the group selection function is the highest, divided 
by the comoving volume occupied by the groups. The 
mock catalogs are used to estimate the effects on the 
observed number density due to our slitmask target se- 
lection (which decreases the number of observed groups, 
as most groups have only two or three observed galax- 
ies) and the false interloper rate due to the group-finder. 
Correcting for these effects, the actual comoving num- 
ber density is estimated to be n = 6 x 10 -4 h 3 Mpc -3 . 
This corresponds to a mean inter-group spacing of d = 
11.8 ft -1 M p c. Fo r comparison to the sample used by 
iGerke et all (|2005), the number density of groups with 
cr > 350 km s -1 is n = 2.4 x 10 -4 h 3 Mpc -3 , after ap- 
plying the above corrections, which corresponds to an 
intergroup spacing of d — 16 hr 1 Mpc. 

For a ACDM cosmology with fl b = 0.05, A = 0.7, 
fi TO = 0.3, h = 0.7, erg = 0.9, a comoving number den- 
sity of n = 6 x 10 -4 /i 3 Mpc -3 results in M min = 5.9 x 
10 12 /i -1 Af Q at z = 0.8 and M m j n = 5.5 x lO 12 /i -1 M 
at z = 1 for a iSheth fc Tormenl 0999) mass function. 
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A comoving density of n = 2.4 x 10 4 h 3 Mpc 3 , es- 
timated for the a > 350 km s" 1 sample, corresponds 
to M mm = 1.2 x I0 13 h- 1 M Q at z = 0.8 and M min = 
1.1 x 10 13 /i _1 AfQ at 2 = 1. These masses are only ap- 
proximate as the group number and volume are both just 
estimates. The errors on the minimum dark matter halo 
masses inferred from the observed number densities are 
likely ~ 50%, given cosmic variance errors and the uncer- 
tainties in the corrections made to the number densities 
from the mock catalogs. 

We check for consistency betwe en the observ e d and 
predict ed clustering of these h alos. IMo fc White! (2002) 
use the lSheth fc Torrrienl l|1999l) model to predict the evo- 
lution in the clustering of dark matter halos and find that 
halos of mass M m i n = 5.5 x 10 12 /i _1 Af Q will have a clus- 
tering amplitude of erg = 1.0 at z — 1. Here a$ is defined 
as the standard deviation of halo count fluctuations in 
a sphere of radius 8 h' 1 Mpc; it can be preferable to 
quoting a scale-length, ro, as it removes the significant 
covariance with 7. eg can be calculated from a power- 
law fit to £(r) using the formula, 

where 

72 

J2(7) = (3- 7 )(4- 7 )(6- 7 )27 (9) 

^Peebles! Il980|> . Using the power-law fits to £(r) for 
groups in the DEEP2 data, we find ag = 1.0. in agr ee- 
ment with the pr e dicted value of lMo fc White! ((2002) . 

IKravtsov etaD (|2004f) find using dark matter simula- 
tions with the same concordance cosmology that a num- 
ber density of n = 6 x 10~ 4 h 3 Mpc~ 3 at z = 1 corre- 
sponds to a minimum mass of Mmin = 5.9 x \Q 12 h~^M Pl . 
This value is comparable to the Sheth fc TormerJ (|1999f) 
value quoted above. They predict a clustering amplitude 
of ro = 5.2 h^ 1 Mpc and 7 = 2.16 for these halos, which 
corresponds to erg = 1.05, in good agreement with our 
observed value of <7g for the DEEP2 groups. 

We verify from the mock catalogs that the actual min- 
imum dark matter halo masses for these group samples 
are similar to those estimated above. For the group sam- 
ple defined as having an estimated a > 200 km s _1 , 
the mass distribution has a rough minimum dark mat- 
ter halo mass of M m i n = 2 — 3 x 10 12 /i -1 M Q and the 
a > 350 km s _1 sample has a rough minimum dark mat- 
ter halo mass of M min = 4 - 5 x 10 12 h -1 M©. The mass 
distributions for the group samples in the mock cata- 
logs do not have a very clearly defined lower-mass cutoff, 
however, and these values are roughly half of the values 
quoted above. Even though the galaxies in the mock cat- 
alogs are randomly drawn from the velocity distribution 
of individual dark matter particles, the fact that only a 
small number of galaxies are observed in a single group 
leads to significant scatter between the estimated a from 
the observed galaxies and the actual a and therefore the 
mass of the underlying dark matter particles; this scatter 
likely accounts for the factor of two discrepancy between 
the halos in the mock catalogs and the predictions of 
iSheth fc TormerJ l|1999lh 

5. CLUSTERING OF GALAXIES IN DIFFERENT 
ENVIRONMENTS 



In this section the clustering properties for galaxies in 
groups relative to the full galaxy population and to field 
galaxies are investigated. Unlike for the group sample, 
which has larger Poisson errors and can only be mea- 
sured on scales r > 3 hr 1 Mpc, where redshift-space 
distortions are small, here we are interested in measur- 
ing clustering properties of galaxies on small scales where 
redshift-space distortions are much more significant. In- 
stead of inferring £(r) in real space from measurements 
of £(s) in redshift space, we measure the projected two- 
point correlation function, w p (r p ), from which £(r) can 
be more directly inferred. Color information additionally 
allows us to divide the sample of group galaxies into red 
and blue populations and measure w p (r p ) for each. We 
test the effects of our slitmask target selection algorithm 
and group- finder on these results using mock catalogs. 
Finally, we are able to empirically separate the observed 
w p (r p ) for galaxies in groups into a 'one-halo' and a 'two- 
halo' term, by keeping pair counts where both galaxies 
are in the same or in different groups. 

5.1. Clustering of Full Galaxy Sample and Galaxies in 
Groups in DEEP2 data 

The full galaxy sample here refers to all DEEP2 galax- 
ies between 0.7 < z < 1.0 in the same six pointings 
for which we have group catalogs. The group and field 
galaxy samples are identified using the a > 200 km s _1 
group catalog, and these samples combined make up the 
full galaxy sample. Fig. [5] shows the spatial distribution 
of galaxies in three DEEP2 pointings, with different sym- 
bols for group (open triangles) and field (crosses) galax- 
ies. In the DEEP2 data 39 +/-4% of all galaxies are 
identified as belonging in recovered groups which have 
a > 200 km s _1 in the redshift range 0.7 < z < 1.0, 
where the error quoted is the standard error across the 
six pointings. This rate is artificially increased by false 
interlopers, but it is also decreased by our slitmask target 
selection by roughly the same amount, as estimated using 
mock galaxy catalogs. In the mock catalogs, 27% of all 
galaxies are identified as belonging in recovered groups 
after target selection (24% are in real groups before tar- 
get selection), significantly less than in the DEEP2 data. 
This difference is most likely due to the free parameters 
in the group-finding algorithm having been tuned to re- 
produce the observed n(a, z) for a > 350 km s" 1 , not the 
a > 200 km s" 1 cutoff in the present sample (see Fig. 
6 in iGerke et alJ ((2005) for details). The clustering re- 
sults shown here do not depend on the absolute number 
of galaxies in each of the samples (all, group, and field 
galaxies) . 

Fig. ^presents w p (r p ) for the full galaxy sample (solid 
lines), galaxies in groups (dashed lines) and for field 
galaxies (dotted lines), in both the DEEP2 data (top) 
and in the mock catalogs (bottom). The top left panel 
compares w p (r p ) as measured in the observed DEEP2 
data (thin lines with no error bars) with w p (r p ) after cor- 
recting for effects due to target selection and the group- 
finder (thick lines with error bars). To make this cor- 
rection, we have used the ratio of w p (r p ) as a function 
of scale in the mock catalogs for group galaxies, field 
galaxies, and all galaxies separately, identified using real 
groups before target selection (thick lines with error bars 
in the bottom right panel), to w p (r p ) for galaxies iden- 
tified using recovered groups after target selection (thick 



lines with error bars in the bottom left panel) . Power-law 
fits to the corrected w p (r p ) points are shown in the top 
right panel of Fig. Eland are listed in tabled along with 
fits to the observed points. Field galaxies are well fit by 
a power-law on scales r p = f — 20 hr 1 Mpc. On smaller 
scales the correlation function is negative, as field galax- 
ies are not found within ~ 0.6 h^ 1 Mpc of each other. 
Pairs of galaxies within that distance are likely to be part 
of a group. 

We will present updated results for w p (r p ) for galax- 
ies in the DEEP2 dataset as a function of galaxy color, 
luminosity, redshift, etc. in a future paper (Coil et al. 
2005, in preparation). Here we focus on the difference 
between the clustering of galaxies in groups relative to 
the full galaxy sample (throughout this paper we use the 
terms "all galaxies" and "the full galaxy sample" inter- 
changeably; they include both field and group galaxies). 
The full galaxy sample used here includes a total of 9787 
galaxies in the redshift range 0.7 < z < 1.0 in the same 
six pointings used for the group analysis. This repre- 
sents a great advance over the sample used in lCoil et al.l 
l)2004aj) . which contained 2219 galaxies in one pointing 
over the redshift range 0.7 < z < 1.35. The fits for tq 
and 7 here agree with our earlier results, but the errors, 
both Poisson and cosmic variance, are much smaller in 
the current sample. For example, we can now address 
whether the signifi cant rise in slope of w n (r n ) on small 
scales predicted bv lKravtsov et all ((2004) for galaxies at 
z = 1 based upon their simulations is present in the 
DEEP2 data. They find that for a galaxy sample with a 
comoving number density of n — 1.5 x 10~ 2 h 3 Mpc~ 3 
(simila r to the DEEP2 sample at z ~ 0.8; see iLin et al.l 
(2004)) that the slope of £(r) changes from 7 ~ 1.65 
when measured on scales of r = 0.3 — 10 h" 1 Mpc, where 
ro is measured to be ~ 4 /i -1 Mpc, to 7 ~ 1.9 over scales 
r ~ 0.1 - 10 hr 1 Mpc, where r ~ 3.5 hr 1 Mpc. Fit- 
ting our corrected w p {r p ) for the full galaxy sample over 
the same range in r p results in r p — 0.3 — 10 hr 1 Mpc, 
r = 3.64 ± 0.07 h' 1 Mpc and 7 = 1.73 ± 0.04, while 
for scales r p = 0.01 - 10 hr 1 Mpc, r = 3.64 ± 0.07 hr 1 
Mpc and 7 = 1.73 ± 0.03, with no change at all in ei- 
ther amplitude or slope. Therefore no evidence is found 
for a rise in the slop e on small scales as predicted by 
iKravtsov et alJ l)2004|) . We also do not find a significant 
difference in the slope on scales r p < 1 h^ 1 Mpc com- 
pared to scales r p > 1 h~ l Mpc, as is seen for the galaxies 
in groups. Fitting for a power-law on scales r p — 0.05 — 1 
h' 1 Mpc, r = 3.52 ±0.16 h^ 1 Mpc and 7 = 1.78 ±0.06, 
while on scales r p = 1 — 20 h^ 1 Mpc, ro = 3.70±0.10 and 
7 = 1.77±0.06. The full galaxy sample therefore appears 
to be consistent with a single power-law slope over the 
range r p = 0.1 — 20 hr 1 Mpc. We note that in the mock 
catalogs, the slope measured for w p (r p ) is consistent with 
the data, though the amplitude of ro is ~ 7% higher. As 
stated before, the mock catalogs were designed to match 
the previously published measurements of w p (r p ) for the 
full DEEP2 galaxy sample. 

The correlation function for galaxies in groups is rel- 
atively well-fit by a power-law over all scales; a broken 
power-law fit with a break at r p = 1 h^ 1 Mpc results 
in a steeper slope on small scales, with low significance. 
The slope for a single power-law is 7 = 2.12 ± 0.06 for 
scales r p = 0.05 — 20 h^ 1 Mpc, while it increases to 



7 = 2.16 ± 0.11 for scales r p = 0.05 - 1 h' 1 Mpc and 
decreases to 7 = 2.02 ± 0.15 for scales r p = 1 — 20 hr 1 
Mpc. Note the significantly different shape of w p (r p ) 
for group galaxies in the mock catalogs, which exhibit a 
strong break at r p ~ 1 h^ 1 Mpc. This sharp rise on small 
scales is not seen in the DEEP2 data. As we will show 
in the next subsection, this difference between the mock 
catalogs and the data is not due to any systematic effect 
from our slitmask algorithm or in our group identifica- 
tion, as the general shape of the correlation function in 
the mock catalogs is not changed by these. This differ- 
ence in the clustering of galaxies in groups between the 
mock catalogs and the data is the first indication that 
the mock catalogs, which were constructed to match the 
z ~ 1 luminosity function and clustering of all galaxies in 
the DEEP2 data, do not reproduce additional properties 
of the data. It appears that the spatial distribution of 
galaxies in groups is less concentrated in the data than 
in the mock catalogs. We discuss the implications of this 
in §9. 

We also investigate wJr p ) for group galaxies as a 
function of color. Fig. shows w p (r p ) for galaxies in 
the DEEP2 data which are identified as belonging to 
groups and have red or blue colors, defined by the min- 
imum in the color bi- modality in restframe (U — B)q, 
at (U — B)q = 1.05. In the group galaxy sample, 20% 
of group galaxies are red, while for the full galaxy sam- 
ple 15% of galaxies are red. Corrections for target se- 
lection and our group-finder are made using the mock 
catalogs as above. Red galaxies in groups have a steeper 
slope in w p (r p ) and a higher correlation length than blue 
galaxies; power-law fits result in ro = 4.77 ± 0.20 ft. -1 
Mpc and 7 = 2.15 ± 0.05 for blue group galaxies and 
r = 5.81 ± 0.45 hr 1 Mpc and 7 = 2.27 ± 0.11 for red 
group galaxies. The steeper slope for the red galaxy sam- 
ple implies that red galaxies are more centrally concen- 
trated in their parent dark matter halos; we investigate 
this more directly in §6. Colors are not currently in- 
cluded in the mock catalogs and so we can not present 
this measurement for the mock catalogs. 

5.2. Effects of Slitmask Target Selection and 
Group-finder 

As for the clustering of groups in the DEEP2 data, 
mock catalogs are used to quantify systematic errors 
in our measurements of w p (r p ) for each galaxy sample, 
where real groups are used to define galaxies in groups 
or in the field. The bottom left panel of Fig. shows 
w p (r p ) measured for all galaxies, group galaxies, and field 
galaxies in the mock catalogs, before (thin lines without 
error bars) and after (thick lines with error bars) tar- 
get selection. Separate power-law least-squares fits to 
w p {r p ) in the mock catalogs are performed for the full 
galaxy sample and for galaxies in groups and in the field; 
the results before and after target selection are shown in 
Table |21 Field galaxies are only affected by the target 
selection on the very smallest scales, r p < 0.2 h~ x Mpc, 
and there is no significant change to the correlation func- 
tion of field galaxies on the scales over which we measure 
a power-law, r p = 1 — 20 h~ x Mpc. For the full galaxy 
sample, our target selection algorithm causes w p (r p ) to 
be slightly underestimated on small scales (r p < 2 h^ 1 
Mpc) due to our inability to target all close neighbors. 
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Before target selection is applied, a power-law fits well 
for the full galaxy sample on scales r p = 1 — 20 h^ 1 Mpc, 
with a steeper slope on small scales, r p = 0.1 — 1 h^ 1 
Mpc. This change in slope on small scales in the mock 
catalogs is more significant (2.5a vs 1.3c) for the sam- 
ple before target selection than after, due to the smaller 
error bars for the larger sample. We note that the scale 
at which the slope ch anges is larger than predicted by 
iKravtsov et alJ (|2004fl in their simulations, and that this 
difference in slope is not seen in the DEEP2 data, as 
discussed above. 

For galaxies in groups, however, the effects of target se- 
lection are more complicated; it increases their observed 
clustering on all scales. In further tests, we find that tar- 
get selection does not affect the measured clustering of 
galaxies in mock catalogs known to belong in real groups, 
so long as the group membership is identified before tar- 
get selection. However, galaxies can be identified as be- 
longing to a group after target selection only if two or 
more observed galaxies are in that group. We can iden- 
tify only half of the groups after target selection that 
were detectable before target selection, and the groups 
which are preferentially lost are those with a low richness. 
This causes w p (r p ) for galaxies in groups to be greater 
when groups are identified after target selection, as only 
galaxies in the richest, and presumably highest mass and 
most clustered, groups will be included in the observed 
sample. 

The mock catalogs are also used to investigate the ef- 
fect of our group-finder on the measured clustering of 
galaxies in groups, as we want to be sure that our group- 
finding algorithm is not imposing (even indirectly) a pre- 
ferred scale for groups, which could potentially artificially 
cause a change in slope at small scales in w p (r p ) for group 
galaxies. This is tested by comparing the clustering of 
galaxies in real and recovered groups in the mock cat- 
alogs, where both samples have had the DEEP2 target 
selection algorithm applied. The results are shown in the 
bottom right panel of Fig. H3 where both real and recov- 
ered groups are seen to have an inflection in w p (r p ) such 
that the slope rises on small scales (r p < 1 h^ 1 Mpc). 
the group-finder, then, is not imposing this scale on the 
clustering results. We also find that the group-finder ef- 
fectively cancels much of the effect of our slitmask target 
selection, in that w p (r p ) for galaxies in recovered groups 
after target selection is similar to w p (r p ) for galaxies in 
real groups before target selection. The overall correc- 
tion applied to the observed w p (r p ) for group galaxies in 
the DEEP2 data is therefore small. 

Throughout the paper we compare w p (r p ) for the 
DEEP2 data with mock catalogs by applying corrections 
to the observed correlation function to account for ef- 
fects of our slitmask target selection and group-finder 
and compare to results in the mock catalogs before tar- 
get selection for real groups. None of our results change 
if instead we compare results for the observed w p (r p ) in 
the data to results in the mock catalogs before target 
selection for real groups. 

5.3. One- and Two-Halo Terms of the Group Galaxy 
Correlation Function 

As it is known which group each of the group galaxies 
belongs to, we can empirically measure the contribution 
to w p (r p ) from pairs of galaxies in the same or in different 



groups. This is akin to determining the 'one-halo' and 
'two-halo' terms of w p (r p ) in the halo model language, 
where each group is identified with a single dark matter 
halo. The total correlation function is then the sum of 
these two terms: 

£(r) = [l + &/,(»•)] +6 h (r). (10) 

We calculate the 'one-halo' and 'two-halo' correlation 
functions using the following estimators: 
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where DD\ and DD2 are pair counts of galaxies within 
the same group and in different groups, respectively. We 
then sum these along the line of sight (in the it direc- 
tion, to iT ma x) to obtain w Pj ih and w p ^h- The projected 
correlation functions sum such that 

)■ (13) 

Fig. |S| shows w Pt ih{r p ) and w Pt 2h{i~ P ) for galaxies in 
groups in both the DEEP2 data (left) and for real groups 
in the mock catalogs (right). The data have been cor- 
rected for effects due to our target slitmask and group- 
finder algorithms. 

In the DEEP2 data the scale at which the 'one-halo' 
and 'two-halo' terms intersect is r p = 1.0 h^ 1 Mpc; the 
scale in the mock catalogs is r p = 0.5 h~ l Mpc. Exactly 
where the break occurs between the 'one-halo' and 'two- 
halo' terms will presumably depend on the type of groups 
we are probing; larger groups may have this break at a 
larger radius. The change in slope seen in w p (r p ) for 
group galaxies in the mock catalogs is easily understood 
as the scale at which the 'one-halo' and 'two-halo' terms 
equally contribute to w p (r p ). However, the 'one-halo' 
term has a very different shape in the DEEP2 data than 
in the mock catalog, such that galaxies in the data which 
belong to the same group are not as clustered on small 
scales as in the mock catalogs. This is likely due the 
mock catalogs having an incorrect spatial distribution 
of galaxies within dark matter halos on small scales; we 
d iscuss this further i n §9. 

lYang et alJ (2005c) measure w p (r p ) for group galaxies 
in 2dF Galaxy Redshift Survey data at z ~ and find 
that it is not well fit by a single power-law; the 'one- 
halo' term is enhanced relative to the 'two-halo' term 
and there is a rise in w p (r p ) on scales r p ~ 1 — 2 h^ 1 
Mpc. The strength of the rise depends on the abun- 
dance and luminosities of the groups; galaxies in more 
luminous (and presumably more massive) groups have a 
larger 'one-halo' term and a stronger rise in the slope of 
w p (r p ) on small scales. It is only for the full galaxy pop- 
ulation (including field galaxies) that they find a single 
power-law fit to w p (r p ). For our group sample at z <~ 1, 
the shape for w-n ir-n) for galaxies in groups is similar to 
what is seen by lYang et aD (l2005cfl for groups with a 



what is seen by I Yang et alJ puuod) tor groups with a 
comparable number density. lYang et al.l l)2005cj) find a 
small rise in the slope of w p (r p ) on scales below r p = 1 
h~ l Mpc but do not quantify this. Our results at z ~ 1 
appear to be similar to their findings at z ~ 0, though 
with larger errors due to our smaller sample size. 
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6. GROUP-GALAXY CROSS-CORRELATION FUNCTION 

In this section we present the cross-correlation func- 
tion between group centers and the full galaxy sample, 
which is sensitive to the radial profile of galaxies in and 
around groups. As with the group and galaxy correlation 
functions, to avoid redshift-space distortions we measure 
the projected cross-correlation, w p (r p ). As discussed in 
§2, errors in the positions of group centers will have some 
effect on scales r p < 1 h~ x Mpc; for this reason we do not 
plot results for r p < 0.3 hr 1 Mpc. However, comparisons 
between data and mock catalogs are unaffected, as the 
mock catalogs have been treated in an identical manner 
as the data. 

Fig. |5] shows the projected cross-correlation between 
group centers and the full galaxy sample in the DEEP2 
data (left) and the mock catalogs (right). The left panel 
shows the observed w p (r p ) as a dashed line and the cor- 
rected w p {r p ) as a solid line, where corrections for target 
selection and the group-finding algorithm as a function 
of scale are made using the ratio of w p (r p ) in the mock 
catalogs between real group centers and all galaxies be- 
fore (solid line, right panel) target selection and between 
recovered group centers and all galaxies after (dashed 
line, right panel) target selection. The dotted line in 
the right panel shows the cross-correlation between real 
groups and all galaxies after target selection. The tar- 
get selection algorithm has the effect of increasing the 
cross-correlation on scales r p > 0.4 hr 1 Mpc, while de- 
creasing the amplitude on smaller scales. The small-scale 
decrease is due to our inability to target galaxies which 
are in close projection on the plane of the sky; this causes 
us to undersample close neighbors. On large scales, the 
effect is due to the slitmask target algorithm affecting 
which groups we identify; after target selection we lose 
many of the pairs of galaxies which were identified as 
groups before such that we preferentially identify the 
groups with more observed members, which are presum- 
ably more massive and therefore more clustered. The 
effect of the group-finding algorithm is to increase the 
cross-correlation on scales r p > 0.4 Mpc, where the 
group-finder has by definition targetted overdensities in 
the galaxy distribution. 

Comparing the solid lines in the two panels of Fig. 
which shows w p (r p ) for real groups before target selec- 
tion to the corrected data, the overall shape of the cross- 
correlation agrees reasonably well, though the amplitude 
is somewhat higher in the mock catalogs on both small 
and large scales. This is consistent with what is seen for 
the correlation function of galaxies in groups shown in 
Fig. which are more strongly clustered in the mock 
catalog than in the DEEP2 data. 

We also investigate the dependence of the radial dis- 
tribution of galaxies in groups on galaxy color. Fig. ITU1 
shows the projected cross-correlation function between 
either red or blue galaxies and group centers, where again 
the galaxy sample has been split at the bi-modality in the 
restframe (U — B) color distribution at (U — B) = 1.05. 
Within groups, on small scales, r p < 0.5 h" 1 Mpc, 
red galaxies are much more strongly clustered than blue 
galaxies, i.e., red galaxies are preferentially found near 
the centers of groups. In the DEEP2 data 20% of group 
galaxies in our sample are red, while 13% of field galaxies 
and 15% of the full galaxy sample (used in this cross- 



correlation) are red. 

Similar trends are seen at z ~ hv lYang et aH<|2005h|h 
who measure the cross-correlation between group centers 
and all galaxies in 2dF and SDSS data. They also find 
a difference between the radial distribution of red and 
blue galaxies, though it is only apparent for groups with 
masses M > 10 13 /i _1 Mq, and the differences are smaller 
than those found here at z ~ 1. 

7. RELATIVE BIAS BETWEEN GROUPS AND GALAXIES 

Measuring the clustering properties of groups, all 
galaxies, and galaxies in groups in the DEEP2 data al- 
lows us to measure the relative bias between galaxies in 
groups and all galaxies and between groups and galaxies. 
Fig. ^2 plots the relative bias of group galaxies to the 
full galaxy sample, which we define as the square root 
of w p (r p ) for group galaxies (dashed lines in Fig. 6) di- 
vided by w p (r p ) for the full galaxy sample (solid lines 
in Fig. 6), as a function of scale for r p = 0.1 — 20 h -1 
Mpc in both the DEEP2 data (top left panel) and the 
the mock catalogs (bottom left panel), after correcting 
for target selection and the group-finder. The bias seen 
between group galaxies and all galaxies is not surprising, 
as group galaxies reside in more massive dark matter ha- 
los. There is a clear scale-dependence to the relative bias 
between group galaxies and all galaxies in the DEEP2 
data, which falls from b re i ~ 2.5 ± 0.3 at r p = 0.1 h^ 1 
Mpc to b rei ~ 1 ± 0.5 at r p = 10 h^ 1 Mpc. The mock 
catalogs have a much higher relative bias on small scales 
( r v 1 Mpc) which does not match the bias seen 
in the data. This reflects the strong rise in slope of the 
correlation function of galaxies in groups seen on small 
scales in the mock catalogs. 

The ratio of the group center-full galaxy sample cross- 
correlation function (Fig. 9) to the galaxy correlation 
function (solid lines in Fig. 6) provides a measure of the 
relative bias of groups to galaxies, which is shown on the 
right side of Fig. ^2 We have corrected for slitmask 
target selection and the group-finder. There is some 
scale-dependence to the relative bias between groups and 
galaxies in the DEEP2 data (top right panel) and the 
weighted mean relative bias is b re i = 1.17 ± 0.04 over 
scales r p = 0.5 — 15 h^ 1 Mpc. The mock catalogs have 
a mean value of b re i = 1.23 ± 0.02 for r p = 0.5 — 15 h" 1 
Mpc, in reasonable agreement with the data, though the 
mock catalogs again show a higher bias on small scales, 
below r p ~ 1 h^ 1 Mpc, and the agreement with the data 
is better on scales r p > 1 h^ 1 Mpc. These measures 
of the relative biases of groups to galaxies and galax- 
ies in groups to all galaxies at z ~ 1 are further con- 
straints which simulations and galaxy evolution models 
must match, in addition to measures of w p (r p ) for all 
galaxies and for groups. We discuss the implications of 
these differences between the data and mock catalogs in 
the next two sections. 

8. EFFECT OF VARYING THE HALO MODEL 
PARAMETERS 

The differences seen between the clustering of group 
galaxies on small scales in the DEEP2 data and the mock 
catalogs could be due to the mock catalogs having the 
wrong spatial distribution for galaxies within their par- 
ent dark matter halos and/or the wrong HOD, which 
specifies the probability that a dark matter halo of mass 
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M hosts N galaxies, P(N\M). To illustrate how much 
these differences may be due to the HOD used to cre- 
ate the mock catalogs, we investigate the clustering of 
group galaxies in two mock catalogs with similar num- 
ber densities and different HODs. In addition to the 
mock catalog used throughout this paper (labeled as 
"B256"), we also analyze a mock catalog in which a dif- 
ferent HOD was applied to the same dark matter simula- 
tion; this model is labeled as "C256" and was chosen as 
one of the most discrepant HOD models that has an ob- 
served w p (r p ) for the full galaxy sample that, by design, 
matches the resu lts for the DEEP2 data published in 
ICoil et alJ (|2004alh The HODs for galaxies with L > L* 
for these two models are shown in the upper left of Fig. 
1121 where model B256 is seen to have a lower minimum 
halo mass hosting a single galaxy, and a steeper slope 
on larger mass scales (~ 0.5 compared to ~ 0.26 for the 
C256 model), which results in having a greater fraction of 
galaxies residing in massive halos. The curves for galax- 
ies with lower luminosity thresholds have a similar shape 
and higher amplitude than wha t is shown here (see Fig. 
1 in lYan. White, fc Coill (|2004)). Mock catalogs made 
with the C256 model have 35% of galaxies in recovered 
groups after target selection, similar to what is found 
for the DEEP2 data (39%), and higher than the value 
found in the B256 model (27%), even though the C256 
model has relatively fewer galaxies in more massive halos. 
This is due to the parameters of the group-finding algo- 
rithm having been tuned to match the observed n(a, z) 
of the DEEP2 data for a > 350 km s _1 ; we have not re- 
tuned the group-finder to the C256 mock catalogs or our 
a > 200 km s^ 1 cutoff. Both the B256 and C256 mock 
catalogs have the same number density for the full galaxy 
sample. The clustering measures shown here have all 
been corrected for slitmask target effects and the group- 
finder. 

Fig. El shows the correlation function for all galax- 
ies (top right panel) and for group galaxies (bottom left 
panel) in each of the two halo model mock catalogs. The 
w p (r p ) for the full galaxy sample is very similar in the 
two catalogs; the only differences are on scales less than 
r p ~ 0.5 hr 1 Mpc, where the C256 model has a slightly 
higher clustering amplitude. For w p (r p ) for group galax- 
ies, the overall shape is similar for the two models but the 
amplitude in model B256 is higher at all scales, as this 
model has more galaxies in massive halos, such that the 
group galaxies will be more clustered. These figures show 
that the amplitude of the group galaxy correlation func- 
tion adds an additional constraint on the HOD, which is 
not gained from the correlation function of all galaxies 
alone. It also shows that the general shape of the group 
galaxy correlation function, and in particular, the rise in 
slope on small scales, is not sensitive to the parameters 
of the HOD used. 

The bottom right panel of Fig. ^| presents the cross- 
correlation function of group centers and all galaxies in 
both mock catalogs. Here there is a difference in the 
shape of the cross-correlation on small scales for the dif- 
ferent HODs. The C256 model has a lower amplitude 
than the B256 model over almost all scales but shows a 
distinct rise on the smallest scales, r p < 0.5 h^ 1 Mpc, 
which is not seen in the B256 model. Indeed, both the 
correlation function for all galaxies and the group-galaxy 
cross-correlation function in the C256 model show a rise 



on small scales that is not seen in the B256 model; this 
results from the C256 model having preferentially more 
galaxies in smaller mass halos which dominate the pair 
counts at small separations. 

Results from the DEEP2 data are also compared to the 
different halo model mock catalogs in Fig. By de- 
sign, w p (r p ) for all galaxies matches both mock catalogs 
well, though the data do not show the rise on the small- 
est scales that is seen for the C256 model. The shape of 
the correlation function for group galaxies in the DEEP2 
data does not match either of the halo model mock cat- 
alogs; the data show a significantly shallower slope on 
small scales. The amplitude of the cross-correlation func- 
tion agrees better with the B256 model than the C256 
model, and the shape of the cross-correlation disagrees 
with the C256 model on small scales. The significant 
difference in w p (r p ) for group galaxies on small scales is 
presumably due to a difference in the spatial distribution 
of galaxies in groups in the data and the mock catalogs, 
as it docs not appear to be reconcilable by altering the 
HOD. This implies that the spatial distribution of galax- 
ies in groups in the mock catalogs is incorrect. This will 
be discussed further in the next section. 

We note that if the C256 mock catalogs had been used 
to correct the observed w p (r p ) for the full DEEP2 galaxy 
sample and galaxies in groups as presented in Fig. [S]for 
slitmask target effects, none of our conclusions in the 
paper would change, as the relative differences before 
and after target selection are similar in the two mock 
catalogs. The differences for all DEEP2 galaxies if using 
the C256 mock catalogs to correct for target selection 
effects are negligable, well within the la errors quoted 
in Table |21 The differences for group galaxies are within 
the 2<t errors, with both the corrected tq and 7 being 
lower (r = 4.72 ± 0.23 h~ l Mpc and 7 = 2.05 ± 0.06), 
and there is still no significant difference in the slope of 
w p {r p ) on small and large scales in the DEEP2 data. 

9. DISCUSSION AND CONCLUSIONS 

Groups bridge the gap between galaxies and clusters 
in both mass and scale, and are also the likely locations 
of galaxy mergers. Measurements of the clustering of 
groups can constrain cosmological parameters, and mea- 
surements of the clustering of galaxies in groups can con- 
strain both halo model parameters and the spatial profile 
of galaxies in their parent dark matter halos. Here we 
present the first results on the clustering of groups and 
galaxies in groups at z ~ 1. We measure four types 
of correlation function statistics in the DEEP2 dataset: 
1) the group correlation function, 2) the galaxy correla- 
tion for the full galaxy sample, 3) the galaxy correlation 
function for galaxies in groups, and 4) the group-galaxy 
cross-correlation function. The first clustering measure 
probes the dark matter halo-halo correlation function on 
mass scales of galaxy groups, which is well-understood 
from dark matter simulations alone. The clustering of 
groups in the DEEP2 data at z ~ 1 matches predictions 
for a ACDM cosmology and is used to estimate the typ- 
ical dark matter masses of the halos the groups studied 
reside in. 

The second clustering measure, the galaxy correlation 
function for the full DEEP2 galaxy sample, is an up- 
date on results u sing early DEEP2 data presented in 
ICoil et ahl l)2004a|) . Here we present this measurement for 
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the full galaxy sample using a much larger dataset, with 
over four times as many galaxies covering three fields 
in the sky; the statistical error on rg is now 2%. We 
show that w p (r p ) for the full galaxy sample provides con- 
straints on the halo occupation distribution (HOD), the 
number of galaxies that reside in a dark matter halo of a 
given mass) , though not in an unique way. There is some 
leeway in how halos can be populated which results in a 
correlation function that matches our measurements, as 
shown in the upper right panel of Fig. 12. 

The third clustering measure, the correlation function 
of galaxies in groups, is similar to the second but is re- 
stricted to galaxies in more massive halos, as measuring 
the clustering of galaxies in groups is sensitive to a higher 
halo mass range than for the full galaxy sample. The con- 
tribution from the 'one-halo' term is necessarily higher 
for galaxies in groups as these galaxies are identified as 
belonging in halos with several other galaxies. This pro- 
vides further constraints on the halo model parameters 
than those from the measurement of £ (r) for all galaxies 
alone. We also find that red galaxies in groups have a 
steeper slope and higher clustering amplitude than blue 
galaxies in groups. 

The fourth clustering measure, group-galaxy cross- 
correlation function, reflects the spatial distribution of 
galaxies within dark matter halos above a given mass, 
and depends as well upon the parameters of the halo 
model. We find using the group center-galaxy cross- 
correlation function that red galaxies are found prefer- 
entially in the centers of groups comp ared to blue galax- 
ies, which has also been seen locally ijCollister fc Lahavl 
120051 lYang et al.ll2005b|) . We find that this trend is in 
place at z ~ 1. 

All four of these measurements are compared to mock 
catalogs constructed from N-body simulations and de- 
pend differently on, and can therefore simultaneously 
constrain, both parameters of the halo model and the 
spatial profile of galaxies within halos. Comparing these 
clustering meas urements in the DEEP2 d ata with the 
mock catalogs of lYan. White, fc ColH l|2004j) . three of the 
four measures agree fairly well with the simulations, with 
the exception of the correlation function for galaxies in 
groups. The clustering amplitude for the full galaxy 
sample roughly matches the mock catalogs; this is by 
design: the catalogs were constructed with an HOD 
that is consistent with earlier DEEP2 clustering results 
for all galaxies. The HOD used is not uniquely deter- 
mined however; the observed w p (r p ) for the full DEEP2 
galaxy sample can be matched with substantially differ- 
ent HODs (two samples are shown in the upper left panel 
of Fig. 12). We leave an improved HOD reconstruc- 
tion from w p (r p ) for the full galaxy sample for a future 
paper, where we will study the clustering as a function 
of galaxy properties such as luminosity, color, redshift, 
etc., in volume-limited samples; here we focus on com- 
paring the clustering of galaxies in groups to all galaxies 
and the different constraints they provide on the HOD. 
We do note that w p (r p ) for all galaxies is fit by a sin- 
gle power-law on scales r p = 0.05 — 20 h^ 1 Mpc, with 
r = 3.63 ± 0.07 and 7 = 1.74 ± 0.03. 

While the mock catalogs have similar projected clus- 
tering for the full galaxy population to the DEEP2 data 
by design, there is a strong discrepancy in the clustering 



of galaxies in groups. The DEEP2 data do not show a 
significant rise in the slope of w p (r p ) on small scales, for 
either group galaxies or the full galaxy sample (upper 
panels of Fig. 6). In contrast, our mock catalogs have 
a very strong rise on small scales for w p (r p ) for group 
galaxies, though not for the full galaxy sample (bottom 
panels of Fig. 6). To test whether this discrepancy can 
be accounted for by the halo model parameters used, we 
analyze mock catalogs constructed with a different HOD 
but similar clustering for the full galaxy sample (Fig. 
12). We find that there is still a rise in slope for the clus- 
tering of galaxies in groups in the second mock catalogs 
(model C256) which is not seen in the data. This result 
is unaffected by our definition of the group center. 

The slope of w p (r p ) for group galaxies on small scales 
should depend quite sensitively on the spatial distribu- 
tion of galaxies within dark matter halos. We therefore 
conclude that there is a difference in the spatial distribu- 
tion of galaxies within their parent dark matter halos in 
the DEEP2 data and our mock catalogs. The mock cat- 
alogs assume no spatial bias, except for the assumption 
of a central galaxy; the most luminous galaxy in a halo 
is placed at the center of the halo, while all subsequent 
galaxies are assigned to random dark matter particles, 
following an NFW profile. Assuming that the brightest 
galaxy occupies the very center of the halo is likely not to 
be correct, as groups at z = 1 are not expected to have 
a large, dominant, bright galaxy in their centers. This 
assumption will result in a higher correlation function on 
small scales. 

In our mock catalogs the satellite galaxies are assumed 
to follow the same NFW profile as the dark matter par- 
ticles; this appears to not be the correct spatial profile 
for the galaxy population at z = 1. There is evidence in 
both simulations and data at z ~ that satellite galax- 
ies do not follow the same spatial profile as the dark 
matter particles. Simulations have found that subha- 
los have a shallower spatial profile than the dark mat- 
ter particles at z ~ (e.g.. iDiemand. Moore, fe Stadel 
l200llC T ao et aUl2004t iNagai fc Kravtsovll2005D . though 
exactly how galaxies are related to subhalos is still not 
entirely known. Observationally, several authors have 
measured the spatial profiles of galaxies in groups in 
data at z ~ 0. Using the Two Micron All Sky Survey, 
iLin. Mohr. fe Stanford! l|2"004|) stack groups and clusters 
to measure the radial mass-to-light profile and find that 
galaxies are less concentrated in the centers of groups 
and c lusters than the dark matter. iCollister fc L ahav 
( 2005) use the 2dF 2PIGG group catalog to directly mea- 
sure the radial profile of galaxies within groups and find 
that galaxies are less centrally concentrated than what 
is seen for dark mat ter particles in simul ations. Similar 
results are fo und bvlHansen et alJ (I2005T) for clusters in 
SDSS, and hv lDiaz et al.l II2005D and lYa.ng et alJ (2005b) 
for groups in SDSS and 2dF. 

The cross-clustering between groups and galaxies 
matches the mock catalogs well on scales r p > 0.5 h" 1 
Mpc; on smaller scales the mock catalogs have a slightly 
steeper slope. The cross-correlation between groups and 
galaxies is linearly proportional to the radial distribu- 
tion of galaxies in groups, while the correlation of group 
galaxies is proportional to the second power of the radial 
distribution; this may be why the shape agreement be- 
tween the data and the mock catalogs is better for the 
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cross-correlation than for the correlation of group galax- 
ies. The group-galaxy cross-correlation function can also 
be affected by uncertainties in the location of the center 
of each group which can dilute the signal on small scales, 
unlik e for t h e corre lation function of galaxies in groups. 
lYang et alJ l)2005bl) find that the cross-correlation be- 
tween group centers and galaxies at z = 0.1 in 2dF and 
SDSS data is lower on scales r p < 0.1 h~ Y Mpc than in 
their mock catalogs, which do not have a spatial bias with 
respect to the dark matter distribution. They create a 
series of mock catalogs with NFW profiles for the galax- 
ies with lower concentration parameters, c, than in the 
dark matter and find that catalogs with concentration 
values of about one-third the value for the dark matter 
halos match their data well. These mock catalogs show 
the same trend which is required here, namely a lower 
cross-correlation on small scales of r p < 0.1 h~ x Mpc. 

We show here that the clustering properties of galaxies 
in groups can be used to break degeneracies among differ- 
ent HODs that can not be distinguished by the clustering 
of all galaxies alone. We find that galaxies in the DEEP2 
data do not have the same spatial profile as in our mock 
catalogs, which assumes a central galaxy surrounded by 
satellite galaxies following an NFW profile. Using the 
clustering statistics presented in this paper will allow us 
to now construct more realistic mock catalogs for the 
DEEP2 survey that have a better constrained HOD and 
radial profile for galaxies within dark matter halos. 
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Table 1. Power law fits to the group correlation functions £(s) and 
£(r) for groups in the deep2 data and in mock catalogs on scales 
r = 3 - 20 h' 1 Mpc. 



Sample 


No. of 


Observed 




Corrected 






groups 


so 


7 


ro 


7 






h- 1 Mpc 




h- 1 Mpc 




DEEP2 data 












N > 2 


460 


7.3 ±0.6 


1.4 ±0.2 


6.2 ±0.4 


1.5 ±0.2 


N > 4 


204 


9.5 ±0.9 


1.6 ±0.3 


8.3 ±0.8 


1.6 ±0.3 


Mock Catalog 












Recovered 


395 


7.9 ±0.5 


1.6 ± 0.2 






Real 


281 


7.5 ±0.8 


1.6 ± 0.3 


6.0 ±0.7 


1.8 ±0.3 


Before Target Selection 


527 






7.4 ±0.4 


1.6 ±0.2 



Table 2. Power law fits to the galaxy correlation function w p (r p ) 

FOR GALAXIES IN THE DEEP2 DATA AND IN MOCK CATALOGS. OBSERVED VALUES 
ARE AFTER TARGET SELECTION, WHILE CORRECTED VALUES ARE BEFORE TARGET 
SELECTION. 



Sample 


No. of r p range Observed 


Corrected 




galaxies ro 7 


r 7 




h^ 1 Mpc h^ 1 Mpc 


h- 1 Mpc 



DEEP2 data 



Group galaxies 


3840 


0.05 


- 20 


5. 


68 


± 


0. 


26 


2 


.12 


± 


0. 


05 


5 


.26 


± 


0. 


27 


2 


12 


± 


0. 


06 






0.05 


- 1 


7. 


21 


± 


0. 


.73 


1 


.96 


± 


0. 


06 


5. 


.08 


± 


0. 


62 


2 


.16 


± 


0. 


11 






1 - 


20 


5. 


.50 


± 


0. 


33 


1 


.95 


± 


0. 


19 


5 


.34 


± 


0. 


30 


2 


.02 


± 


0. 


15 


Red group galaxies 


794 


0.05 


- 20 


6. 


31 


± 


0. 


44 


2 


.29 


± 


0. 


09 


5 


.81 


± 


0. 


45 


2 


.27 


± 


0. 


11 


Blue group galaxies 


3046 


0.05 


- 20 


5. 


10 


± 


0. 


.20 


2 


.12 


± 


0. 


04 


4 


.77 


± 


0. 


20 


2 


.15 


± 


0. 


05 


All galaxies 


9787 


0.05 


- 20 


3. 


46 


± 


0. 


.09 


1 


.68 


± 


0. 


04 


3 


.63 


± 


0. 


07 


1 


.74 


± 


0. 


03 


Field galaxies 


5947 


1 - 


20 


2. 


55 


± 


0. 


26 


1 


.68 


± 


0. 


12 


2 


.76 


± 


0. 


20 


1 


.72 


± 


0. 


12 


Mock Catalogs (B256) 
















































Group galaxies 


3020 


0.05 


- 1 


3. 


63 


± 


0. 


.35 


2 


.84 


± 


0. 


.12 


3 


.02 


± 


0. 


25 


2 


.99 


± 


0. 


13 






1 - 


20 


5. 


57 


± 


0. 


37 


1 


.75 


± 


0. 


16 


5 


.43 


± 


0. 


29 


1 


,84 


± 


0. 


11 


All galaxies 


11262 


0.05 


- 20 


3. 


.75 


± 


0. 


09 


1 


.69 


± 


0. 


03 


3 


.75 


± 


0. 


09 


1 


.69 


± 


0. 


03 






0.05 


- 1 


3. 


.55 


± 


0. 


18 


1 


.76 


± 


0. 


07 


3 


.52 


± 


0. 


14 


1 


.85 


± 


0. 


05 






1 - 


20 


3. 


71 


± 


0. 


15 


1 


.65 


± 


0. 


07 


3 


.85 


± 


0. 


11 


1 


.68 


± 


0. 


05 


Field galaxies 


8242 


1 - 


20 


3. 


01 


± 


0. 


14 


1 


.57 


± 


0. 


06 


3 


.35 


± 


0. 


12 


1 


.63 


± 


0. 


05 



15 
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Fig. 1. — Observed richness of groups in six pointings in the DEEP2 data and in mock catalogs, for both real and recovered groups. The 
distributions roughly decrease as a power-law, with the largest groups having -~10 members. There are significantly more recovered groups 
than real groups in the mock catalogs. 




1.0 
redshift 



Fig. 2. — Redshift histogram of groups with a > 200 km s~ 1 in the DEEP2 data, using our six most complete pointings to date. The 
solid line is a smoothed version of the histogram which is used to estimate our selection function as a function of redshift for groups in our 
data; for the clustering analyses presented here we only use groups and galaxies in the redshift range 0.7 < z < 1.0. 
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Fig. 3. — Spatial distribution of groups with a > 200 km s 1 in three of the six DEEP2 pointings used, as a function of redshift (or 
comoving line of sight distance) for 0.7 < z < 1.0 and transverse distance. We have projected through the short dimension of each field, 
corresponding to ~ 20 h~ 1 Mpc. The aspect ratio is not true in this figure, such that large-scale features appear compressed along the line 
of sight. The groups generally trace the largest structures and are often seen along filaments. 




Fig. 4. — The clustering of groups with a > 200 km s _1 in DEEP2 data and mock catalogs. The redshift-spacc correlation function, 
£(s) (left), and real-space correlation function, £(r) (right), for the groups with four or more observed members (TV > 4, dashed line) are 
higher than for groups with two or more observed members (N > 2, solid line). The average £(s) and £(r) for groups with TV > 2 in the 
mock catalogs (dotted line) is shown as well, and agrees well with the data. Corrections have been made for the DEEP2 slitmask target 
selection algorithm. 
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Fig. 5. — Spatial distribution of DEEP2 galaxies identified as being field galaxies (thin crosses) or groups galaxies (thick stars) for data 
in the same three pointings plotted in Fig. |2] As in Fig. |3 we have projected through the short dimension of each field, corresponding to 
~ 20 h' 1 Mpc. 



18 




10000 



1000 



100 



0.1 



After Torget Selection 
Before Target Selection 



Mock catalog 




10000 



1000 



100 



Recovered Groups 
Real Groups 



"a. 



Mock catalog 




1.0 

, (h" 1 Mpc) 



10.0 



0.1 



1.0 

, (h" 1 Mpc) 



10.0 



Fig. 6. — Projected two-point correlation function, w p (r p ), for the full galaxy sample (solid lines) and for group galaxy (dashed lines) and 
field galaxy (dotted lines) samples in both the DEEP2 data (top) and mock catalogs (bottom). Top left: Observed correlation functions 
are shown as thin lines without error bars; thick lines with error bars include corrections for the slitmask target selection algorithm and, 
for galaxies in groups or in the field, corrections for the group-finder. Top right: Power-law fits to the corrected Wp(r p ) are shown as 
dot-dashed lines and are given in Table 2. Bottom left: The effects of our slitmask target selection are shown for galaxies in the mock 
catalogs, before (thin lines without error bars) and after (thick lines with error bars) applying the target selection algorithm, where we 
have used real groups to identify group and field galaxies. Bottom right: The effects of the group-finder in the mock catalogs are shown 
for galaxies identified as belonging to real or recovered groups and for field galaxies (there is no difference for the full galaxy sample). 
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Fig. 7. — Projected two-point correlation function, w p (r p ), for DEEP2 galaxies in groups, where we have split the galaxy sample by color 
using the observed bi-modality in restframe (U — B)o- Corrections have been applied for effects due to our slitmask target algorithm and 
group-finder, using the ratio of w p (r p ) in the mock catalogs for galaxies identified as belonging in real groups before target selection (thick 
dashed line in bottom left panel of Fig. 6) compared to galaxies in recovered groups after target selection (thin dashed line in bottom right 
panel in Fig. 6). Red group galaxies (dashed line) as more strongly clustered than blue group galaxies (solid line) and show a steeper 
clustering slope. Power-law fits are plotted as thin dash-dot lines and given in Table 2. 




Fig. 8. — Projected two-point correlation function, w p (r p ), for galaxies in groups in the DEEP2 data (left) and real groups in the mock 
catalogs (right), where the solid and dashed lines include pairs inside the same group ('one-halo') or between groups ('two-halo'). Left: 
Corrections have been applied for our slitmask target selection algorithm and group-finder. The scale at which the two terms intersect is 
somewhat larger in the data than in the mock catalogs, where the amplitude of the 'one-halo' term is significantly higher. 
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Fig. 9. — Projected cross-correlation function, w p (r p ), between group centers and all galaxies in the DEEP2 data (left) and in mock 
catalogs (right). In the left panel, the dashed line shows the observed cross-correlation function, while the solid line shows the function after 
correcting for effects due to target slitmask and the group-finder. In the right panel, the solid and dashed lines show the cross-correlation 
function for real groups before slitmask target selection is applied and recovered groups after the target selection is applied, while the 
dotted line shows the cross-correlation for real groups after target selection. The shape of the cross-correlation in the data agrees well with 
the mock catalogs, though the amplitude is lower. 
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Fig. 10. — Projected cross-correlation function, w p (r p ), between group centers and the full galaxy sample in the DEEP2 data, for galaxies 
rcdward (solid line) and blucward (dashed line) of the observed bi-modality in restframe (U — B)o color. Corrections have been applied for 
our slitmask target selection algorithm and group- finder. Red galaxies are found preferentially near the centers of groups at z ~ 1. 
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r p (h- 1 Mpc) r p (h~ 1 Mpc) 

Fig. 11. — Scale-dependence of the relative bias between galaxies in groups and the full galaxy sample (left) and between groups and 
galaxies (right) in the DEEP2 data (top) and mock catalogs (bottom). The data do not show nearly as strong of a scale-dependence in the 
relative bias between group galaxies and all galaxies on small scales (upper left) as is seen in the mock catalogs (bottom left). The relative 
bias between groups and galaxies in the data (upper right) agrees reasonably well the mock catalogs (bottom right). 
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Fig. 12. — The effects of varying the halo occupation distribution (HOD) parameters on various clustering measures. Top left: HODs 
for the two mock catalogs analyzed here. Shown is the average number of galaxies with L > L* placed in dark matter halos, as a function 
of the halo mass. The projected correlation function for all galaxies (top right) is similar in both of the mock catalogs (dotted and dashed 
lines), which agree fairly well with the DEEP2 data (solid line). The projected correlation function for galaxies in groups (bottom left) 
shows a rise on small scales in both mock catalogs that is not seen in the DEEP2 data. Varying the HOD does not appear to resolve this 
discrepancy. The projected cross-correlation function between group centers and all galaxies (bottom right) agrees fairly well between the 
mock catalogs and data. 



